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ABSTRACT 



Aims. We perform a comprehensive multiwavelength analysis of a sample of 20 starburst galaxies that show a substantial 
»-™H, population of Wolf-Rayet (WR) stars. In this paper, the third of the series, we present the analysis of the O and WR 

star populations within these galaxies. 

Methods. We study the spatial localization of the WR-rich clusters via the detection of the blue WR bump, mainly 
composed by the broad Hen A4686 line and attributed to WN stars, and the red WR bump, composed by the broad 
Civ A5808 line observed in WC stars. We perform a detailed fitting of the nebular and broad emission lines within 
these broad features and derive the numbers of WN, WC and O stars using (i) the standard assumption of constant 
WR luminosities and (ii) considering metallicity-dependent WR luminosities. We then compare our results with the 
predictions given by evolutionary synthesis models and with previous empirical results. 

Results. We report the detection of blue WR bump in 20 regions, but the red WR bump is only detected in six. Aperture 
effects and the exact positioning of the slit onto the WR-rich bursts play a fundamental role in their detection. The 
nebular Hen A4686 line is detected in 21 regions; its intensity clearly decreases with increasing metallicity. We derive 
an empirical estimation of the WNL/(WNL+0) ratio using the intensity of the broad Hen line assuming metallicity- 
dependent WR luminosities. As expected, the total number of WR stars increases with increasing metallicity, but objects 
with 12+log(0/H)<8.2 show a rather constant WR/(WR+0) ratio. The computed WCE/WNL ratios are different 
than those empirically found in nearby star-forming galaxies, indicating that the observed galaxies are experiencing a 
strong and very short burst. Considering metallicity-dependent WR luminosities, our data agree with a Salpeter-like 
IMF in all regimes. 

Conclusions. We consider that the contribution of the WCE stars is not negligible at low metallicities, but deeper 
observations are needed to detect the red WR bump because of the decreasing of the WR luminosities with decreasing 
metallicity. Although available models reproduce the WR properties at high metallicities fairly well, new evolutionary 
synthesis models for young starbursts including all involved parameters (age, metallicity, star-formation history, IMF 
and WR stars properties such as metallicity-dependent WR luminosities, stellar rotation and the WR binnary channel) 
are absolutely needed to perform an appropriate comparison with the observational data. 

O ' , 
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1. Introduction material during their whole lifetime (jMacdcr 198l|) but also 

^ in supernova explosions; finally, they deposit mechanical en- 

' Despite their relatively low number and short lifetime in e rgy, first via strong stellar winds and later as supernovae 

r\ . terms of evolutionary timescales, massive stars have a fun- ([Abbottl fl982T). which are a triggering mechanism of star 

. damental influence on the interstellar medium (ISM) and fo rm ation lWoodwardl[T978l) . Massive stars are the progen- 

galaxy evolution: they generate most of the ultraviolet ion- itors of the most energetic phenomena known today, the 

izing radiation in galaxies, powering the far-infrared lumi- long gamma _ ray bursts (GRB), as they collapse after super - 

nosities through the heating of dust; they enrich the in- nova ex plotions into black holes (1 Wooslev fe Hegeri I2TJ0II . 

terstellar medium, not only by returning nuclear processed The most mass i ve (M > 25 M for Z Q ) luminous (10 5 to 

10 6 Lq) and hot (temperatures ~50 000 K) O stars evolve 

Send offprint requests to: Angel R. Lopez-Sanchez, e-mail: to the Wolf-Rayet phase between 2 and 5 Myr since their 

Angel.Lopez-Sanchez@csiro.au birth. A WR star is interpreted as a central He-burning 

* Based on observations made with NOT (Nordic Optical core that has lost the main part of its H-rich envelope via 

Telescope), INT (Isaac Newton Telescope) and WHT (William strong winds (|Contilll976l lMaederlll990l 119911) . This stage 

Herschel Telescope) operated on the island of La Palma jointly f the most massive stars lasts only some few hundreds of 

by Denmark, Finland, Iceland, Norway and Sweden (NOT) thousands of years (t W R < 5 X 10 5 yr) before they e xplode 

or the Isaac Newton Group (INT WHT) in the Spanish ag T Ib/Ic super novae dMevnet fc Maeded [20051) . 
Observatorio del Roque de Los Muchachos or the Instituto de 
AstroHsica de Canarias. 
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The broad emission features that characterize the spec- 
tra of WR stars are often observed in extragalactic H II re- 
gions. Actually, these Wolf-Rayet galaxies make up a very 
inhomogeneous class of star-forming objects, all having in 
common ongoing or recent star formation which has pro- 
duced stars massive enough to e volve to the WR stage 
(jSchaerer. Contini fe Pindaolll999l ). The blend of the broad 
Hen A4686, Cm/Civ A4650 and Nm A4640 emission lines 
constitutes the blue WR bump; it mainly originates in 
WN stars with a minor contribution of WC stars. The orig- 
inal blue WR bump detectio n was made in the blue com - 
pact d warf galaxy He 2-10 bvlAllen. Wright fe Gossl (fl976t). 
Later, iKunth fe Sargend (|l98lh recognized the importance 
of this WR bump to characterize the mode of star formation 
in starbursts (burst versus continuous mode) and developed 
the first method to empirically derive the WR/O ratio. The 
blend of the Cm A5698 and Civ A5808 broad emission 
lines constitutes the red WR bump. C iv A5808 is the 
strongest emission line in WC stars, but it is barely seen 
in WN stars. The first detections of the r ed WR bump in 
integr ated galaxy spectra were reported bv lKunth fe Schildl 
(|1986D and lDinerstein fe Shields! ()1986D . but it seems to be 
more difficult to detect than the blue WR bump and it 
is always weaker than the blue WR bump ()Guseva et al.l 
I2000L iFernandes et al.ll2004l) . Actually, the red WR bump 
is rarely observed at low metallicities. However, both WN 
and WC stars have been detected in the lowest meta llicity 
{Z ~l/50 Zg) galaxy known, I Zw 18 (llzotov et al.lll997t 
iLegrand et al.ill997t IPeMello et al.|[l998f) . 

Although the main feature of WR galaxies is a broad 
Hen A4686 emission line, a considerable fraction of ob- 
jects also show the nebular (narrow) Hen A4686 line. This 
emission line is very rarely found in Ga lactic H n regions 
(|Garnett et al.lll99ll ISchaerer et al.lll997l ) but common m 
planetary nebulae, and has been sometimes observed in gi- 
ant H II regions where no WR stars have been detected. 
Sources with sufficient photons of energy >54 eV are needed 
to produce the nebular He n emiss i on line, so its origin 
is still puzzling (|Garnett et al.lfl99lUGarnettj|2004l) . Some 
observations concluded that this nebular li ne is intimately 
linke d with the appearan ce of hot WR stars (ISchaerer et all 
[l996t ISchaerer fe Vaccal[l99l iCrowther fe Hadfieldl I200CT 
but recent analyses suggest that O stars may also contribute 
to the He n ionizing flux at low metallicities (Brichmann, 
Kunth & Durret 2008, hereafter BKD08). In one way or 
in another, studies of the nebular Hen emission line give 
clues about the physical processes in H n regions, the ion- 
izing fluxes of starbursts and their contribution to the ion- 
ization of the intergalactic medium (Garnett et al. 1991; 
Scha erer et al. 1998; Sta sihska 1998). This is the reason 
why ISchaerer et al.l (|1999l ) also compiled all the extragalac- 
tic objects (54) showing the nebular Hen A4686 emission 
lin e in their catalog of WR galaxies. The analysis performed 
bv lThuan fe Izotovl (|2005l) listed 465 Blue Compact Dwarf 
galaxies (BCDGs) showing nebular He n A4686 emission. 

Wolf-Rayet stars can be individually observed in our 
Galaxy, the Magellanic Clouds, in some galaxies of the 
Local Group (M 31, M 33, IC 10, NGC 6822, IC 1613) and 
in the brightest galaxies of the Sculptor group (NGC 300 
and NGC 55). A recent re view of the physic al properties of 
WR stars was presented bv lCrowtherl <|2007h . Thanks to the 
analysis of these individual WR stars a quantitative clas- 
sification in WN, WC and WO stars could be established. 
This led the developing of stellar population synthesis mod- 



els tailored to the analysis of massive star populations in 
young starburst galaxies. The best example of these models 
are those presented by Schaerer & Vacca (1998; hereafter 
SV98), which provide detailed predictions of many stellar 
and neb ular UV/optical features. Famous Starburst 99 
models (|Leitherer et al.l 119991: IVazquez fc Leithererl l2005h 
also include predictions for the WR and O stars popu- 
lations. Except for some very low metallicity objects, a 
good agreement between the observations and the synthe- 
sis models of SV98 is usually found. For objects with sub- 
solar metallicity (mostly BCDGs), this comparison indi- 
cates fairly short timescales of star formation, with burst 
ages lower than 2-4 Myr, which is the initial mass function 
(IMF) compatible with a Salpeter slope, but requires the 
existence of high-mass stars (Guseva et al. 2000; De Mello 
et al. 1998; Mas-Hesse & Kunth 1999, Fernandes et al. 2004; 
Buckal ew, Kobulnicky fe Dufour 2005). However, some au- 
thors (|Zhang et al.l I2007D suggested a dependence of the 
slope of the IMF on metallicity. The timescales of star for- 
mation in high-metallicity enviroments are more extended 
than those fo und at low metallicity, with burst duration s 
of 4-10 Myr (jSchaerer et al.l l2000t IFernandes et al"1l2004D . 
although observations are also explained assuming a super- 
position of several bursts. 

We performed a detailed photometric and spectroscopic 
analysis of a sample of 20 WR galaxies. Our main aim 
is the study of the formation of massive stars in star- 
burst galaxies and the role that the interactions with 
or between dwarf galaxies and/or low surface brightness 
objects have in its triggering mechanism. In Paper I 
(|L6pez- Sanchez fc Estebanl 120081 ) we explained the moti- 
vation of this work, compiled the list of the analyzed WR 
galaxies (Table 1 of Paper I) and presented the results of 
optical/ NIR broad-band and Ha photometry. In Paper II 
(Lopez-Sanchez & Esteban 2009) we presented the results 
of the analysis of intermediate-resolution long-slit spec- 
troscopy of 16 objects of our sample of WR galaxies -the re- 
sults for the other four objects were published separately. In 
many cases, two or more slit positions were used to analyze 
the most interesting zones, knots or morphological struc- 
tures belonging to each galaxy or even surrounding objects. 
In this paper we analyze the localization and the proper- 
ties of the O and WR populations within those galaxies. 
Papers IV and V will present the global multiwavelength 
analysis of our WR galaxy sample considering all available 
data. 

In Sect. 2 we describe the detection of the WR fea- 
tures. The procedure used to fit the broad and nebular 
emission lines in the blue and red WR bumps is explained in 
Sect. 3. We analyze the detection of the nebular He n A4686 
line in Sect. 4. The quantitative analysis of the WNL and 
WCE star populations in described in Sect. 5. In Sect. 6 we 
analyze the metallicity dep endence of the WR/(WR +Q) 
ratio. The com parison with ISchaerer fc Vacca] (Il998l) and 
Star burst 99 (jLeitherer et al.lll999UVazauez fc Leithererl 
I2005D models is shown in Sect. 7. We analyze the derived 
WCE/WNL ratio in Sect. 8, while the dependence of the 
IMF on the WR/(WR+0) ratio is discussed in Sect. 9. 
Finally, we compile our main conclusions in Sect. 10. 



2. Detection of the WR features 

As our galaxy sample was extracted from the 
catalog of Wolf-Rayet galaxies compiled by 



2 



Lopez-Sanchez & Esteban: Massive star formation in Wolf-Rayet galaxies III: O and WR populations 
Table 1. Detection of WR features and the nebular Hen emission line for the galaxies analyzed in this work. 



Galaxy 


Zone 


H/3 






He II A4686 


Blue WR bump c 


Red WR bump d 


/ 

(a) 


-W 


///(H/3) 
(b) 


-W 

[A] 


1/1(110) 
00 


-W 
[A] 


///(H/3) 
(b) 


-W 

[A] 


HCG 31 


AC e 


46.13+0.18 


91.1+2.1 


0.25 


n n_i_n q 


U.DDztU.zz 


Q 04-1 Q 


o c:4-n g. 


O A -LCI a 


o. litu.o 




Fl 


3.187+0.016 


218+13 


0.80 


o q I n o 
z.o±U.o 


A Q_l_n A 












F2 


2.28+0.19 


256+43 


0.80 


o o_i_n 


o.y±1.0 


u.yo±u. oo 


*3 A 4-1 P. 






Mkn 1087 


Cf 


78.1 




0.25 






10.6 








Haro 15 


C 


23.3+1.1 


16.4+1.1 


0.40 


0.49: 


0.08: 


2.82: 


0.5: 








A 


23.4+0.9 


75.7+4.2 


0.50 


1.7±0.4 


1.0+0.2 


1.4+0.5 


1.0+0.3 






Mkn 1199 


C e 


74.2+3.1 


21.4+1.3 


0.50 


0.24: 


0.05: 


6.3+1.8 


1.3+0.3 


9.4±2.3 


2.8±0.7 


Mkn 5 


A 


17.7+0.7 


75+5 


0.40 


0.9±0.2 


1.6+0.3 


1.5+0.5 


3.0+1.2 






IRAS 08208+2816 


C 


12.9+0.5 


80+5 


0.25 






4.1+1.3 


3.3+1.0 


3.5±1.3 


3.6±1.4 


IRAS 08339+6517 


#1 


36.5+1.3 


25+2 


0.50 


0.48: 


0.1: 


4.1: 


2.2: 






POX 4 


c 


56.0+1.8 


200+9 


0.60 


1 on i r\ i c 

1.20±0.16 


1.9+0.3 


1.5+0.3 


3.0+0.9 


1.4±0.2 


3.9±0.8 


UM 420 




6.88+0.27 


169+10 


0.50 


1.0+0.2 


2.0+0.4 


1.5+0.5 


2.8+0.7 






SBS 0926+606 


A 


16.2+0.6 


125+6 


0.60 


0.7+0.2 


1.9+0.3 


0.9: 


3.2: 






SBS 0948+532 




8.44+0.32 


213+11 


0.60 


1.3+0.3 


1.6+0.6 


1.1+0.3 


2.1+0.7 






SBS 1054+365 


C 


14.6+0.7 


89+7 


0.40 


0.6+0.2 


0.9+0.3 


1.6+0.5 


2.4+0.8 






SBS 1211+540 




1.84+0.09 


135+10 


0.50 






^ 1 o. 

<l.z. 








SBS 1319+579 


A 


14.6+0.5 


285+14 


0.90 


U.o. 


L.O. 


u. / . 


1 Q. 






SBS 1415+437 


C 


18.5+0.7 


222+11 


0.60 


2.4+0.2 


4.8+0.3 


<0.9: 








III Zw 107 


A 


22.3+0.8 


44.2+2.3 


0.50 






3.1+1.1 


1.3+0.4 






roi y 






oo.o±z.U 


U.oU 


0.33+0.15 


0.24+0.7 


7.9+1.3 


2.5+0.4 


5.8: 


1.9: 


Tol 1457-262 


A 


222+8 


101+6 


0.50 


1.9+0.4 


1.5+0.3 












B 


47.7+2.3 


83+7 


0.40 


2.9+0.3 


2.1+0.2 










ESO 566-7 




10.2+0.8 


13+2 


0.25 














ESO 566-8 




111+5 


95+7 


0.50 


0.74+0.25 


0.65+0.25 


6.2+2.2 


4.1+1.5 


4.8+1.2 


4.2+1.2 


NGC 5253 


A 


134.5+0.4 


234+5 


1.00 


0.10+0.04 


0.21+0.09 


0.71+0.12 


1.6+0.3 


NO OBS. 9 


NO OBS. 9 




B 


135.2+4.2 


254+5 


1.10 


0.05: 


0.11: 


0.71+0.13 


1.7+0.3 


NO OBS. 9 


NO OBS. 9 




C e ' 9 


101.3+3.3 


94.0+1.0 


0.40 


0.06: 


0.05: 


4.13+0.19 


3.55+0.17 


1.00±0.14 h 


NO OBS. 9 




D 


25.6+0.9 


39.0+0.5 


0.30 


0.40+0.11 


0.12+0.04 






NO OBS. 9 


NO OBS. 9 



a In units of 10~ 15 erg s _1 cm -2 and not corrected for extinction. 
6 Assuming /(H/3)=100 and reddening-corrected. 
c Broad He II A4686 emission line only. 
d Broad C IV A5808 emission line only. 

e Broad Nm A4640 also detected. The reddened-corrected ///(H/3) ratios are 0.91, 5.6 and 1.44 for HCG 31 AC, Mkn 1199 and 

NGC 5253 C , respectively. See text for details. 

f Data from lVaceli et all J1997D . see lLopez-Sanchez et al.l i2004bD . 

9 NO OBS indicates that the spectral range where the feature lies was not observed. 

h Broad Civ A4658 detected. The ///(H/3) ratio given here is assuming L wce (A5808)=L WC e (A4658)/1.71 (SV98). See text for details. 



ISchaerer. Contini fc Piiidaol (|1999D . we first analyze 
the presence or absence of the WR features in the star- 
forming bursts within our starburst galaxies and the 
characteristics of their WR populations. Table Q] compiles 
the objects or regions where any WR bump and/or the 
nebular Hen A4686 emission line have been detected. 
Figures 36 and 37 in Paper II show the details of the 
optical spectrum in the 4600-4750 A and the 5550-6000 A 
ranges, respectively, of all important objects (faint regions 
with very low S/N have been excluded). 

The nebular Hen A4686 emission line is detected in 
21 regions within 16 systems (Fig. 36 in Paper II). It is 
unambiguously clear in 14 objects (HCG 31 AC and Fl, 
Haro 15 A, Mkn 5, POX 4, UM 420, SBS 0926+606 A, 
SBS 0948+532, SBS 1415+437 C, Tol 1457-262 A and B, 
ESO 566-8 and NGC 5253 A and D), and particularly 
strong in POX 4. 

The blue WR bump is clearly detected in many objects 
(Fig. 36 in Paper II), and unambiguously observed in 12 re- 
gions (HCG 31 AC, Mkn 1199, Mkn 5, IRAS 08208+2816, 
POX 4, SBS 0948+532, SBS 1054+365, III Zw 107, Tol 9 
and NGC 5253 A,B and C). We indicate the spatial lo- 
calization of the WR stars in each galaxy in Col. 2 of 
Table [I] sometimes they are detected in different regions 



within the same system. Our observations of the galaxy 
IRAS 08339+6517, no cataloged so far as WR, suggest the 
detection of this kind of m assive stars in its central burst 
(jLopez-Sanchez et al.l [2006). Only in three objects previ- 
ously listed as WR galaxies (Mkn 1087, SBS 1211+540 and 
ESO 566-7) do we not detect any feature that can be at- 
tributed to this kind of massive stars. We consider that 
aperture effects and the exact positioning of the slit onto 
the WR-rich bursts play a fundamental role in the detec- 
tion of the WR features. Furthermore, a good S/N ratio is 
needed to observe this feature in faint regions. A good spec- 
tral resolution is also needed to get a proper deblending of 
the nebular and broad He n A4686 lines. 

The dominant contribution to the blue WR bump is 
the broad Hen A4686 line, arising in late type WN stars 
(WNL). Three regions (HCG 31 AC, Mkn 1199 C and knot 
C in NGC 5253) show a clear detection of the broad NlH 
A4640 line, also mainly originated in WNL stars. The ab- 
sence of the faint N v A4604 indicates the predominance of 
WNL stars, although some contribution of early-type WN 
stars (WNE) might be still present in the blue WR bump 
(SV98). 

As was explained in the introduction, the Hen A4686 
emission line and the blue WR bump are mainly linked 
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Fig. 1. Examples of fits to the blue WR bump for some of the analyzed galaxies in the 4620 - 4750 A range. The red dotted 
line represents the position of the He n A4686 emission line, the blue dotted lines indicate the position of [Fe ill] A4658 and [Ar iv] 
AA4711,4740 emission lines. The black dashed line represents the continuum fit. The best fit to the observed spectrum (blue 
continuous line) is shown with a black continuous line, and it is usually composed by the superposition of a broad and a narrow 
Hell A4686 lines. Region D in NGC 5253 only shows a narrow nebular Hell component. The broad Nm A4640 emission feature 
is detected and fitted in HCG 31 AC and NGC 5253 C. This last region also shows a broad ClV A4658 emission line, which was 
included in the fit. The NGC 5253 C also shows the three dashed broad Gaussians that composed the blue WR bump. All diagrams 
include the residual spectrum after subtracting the best-fit model to the observed spectrum. 



to WN stars, while the Civ A5808 emission line (the red 
WR bump) essentially originates in WC stars. As pre- 
vious studies indicate, the red WR bump is much more 
difficult to observe than the blue WR bump. Figure 37 
in Paper II shows the spectral range where the red WR 
bump is located for all important objects; only two galax- 
ies (SBS 1211+540 and NGC 5253) are lacking data in this 
spectral range. There is always a clear detection of the He I 
A5876 line in all spectra, and the faint auroral [N n] A5755 



emission line in also seen in many cases. We clearly identify 
the broad C iv A5808 line in 2 galaxies (HCG 31 AC and 
POX 4), detect it in two galaxies (Mkn 1199 and ESO 566- 
8) and it also seems to be observed in other two galaxies 
(IRAS 08208+2816 and Tol 9). However, we do n ot see this 
feature in galaxies for which iGuseva etaD (|2000h reported 
its detection (HGC 31 G, Mkn 5, UM 420, SBS 0926+606, 
SBS 0948+532, SBS 1054+365 and SBS 1319+579A). The 
non-detection of the red WR bump in these objects may 
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Fig. 2. Examples of fits to the red WR bump for all regions for which we have a detection of the broad Civ A5808 line in the 
5550 - 6000 A range. The red dotted line represents the position of the C IV A5808 emission line, the blue dotted lines indicate 
the position of [Nil] A5755 and He I A5875 emission lines. The black dashed line represents the continuum fit. The best fit to the 
observed spectrum (blue continuous line) is shown with a black continuous line, and it is composed by a broad Civ A5808 and a 
nebular [N n] A5755 line. The broad feature is very evident in HCG 31 AC and in POX 4, detected in Mkn 1199 and ESO 566-8, and 
probably observed (but with a very low S/N ratio) in IRAS 08208+2816 and Tol 9. Except in B£»fe4 and IRAS 08208+2816, the 
spectra also show a clear He I A5875 absorption feature. All diagrams include the residual spectrum after subtracting the best-fit 
model to the observed spectrum. 



be a co nsequence of apertur e effects or the position of the 
slit, i.e.. lGuseva et al.l (|2000t ) observed slightly different re- 
gions than the areas we analyze here. However, our spec- 
tra usually have higher spatial and spectral resolution (see 
Paper II), and even higher signal-to-noise ratio (the av- 
erage SNR value achieved in this part of the spectra for 
our galaxy sample is 10 - 12), particularly those obtained 
using the ISIS spect rograph at th e 4.2m WHT, than the 
spectra compiled bv iGuseva et al.l (|2000D . which were ob- 
tained using 2m and 4m class telescopes (see their Fig. 2). 
In the past, the non-detection of the red WR bump may 
have been a detection-threshold effect (the slope in star- 
bursts rises towards the blue and hence on average there 
are half as many counts in the red WR bump than in the 
blue WR bump), but we consider this not to be the situ- 
ation here because of both the quality of our spectra and 
their relatively good SNR. 

The non-detection of the red WR bump could be ex- 
plained because WC stars are very rarely formed in low- 
metallicity environments, as evolu tionary synthesi s models 
predict and observations suggest (|Crowtherll2007l) . indeed 
our galaxy sample is mainly composed by low-metallicity 
objects. However, in some objects the non-detection of the 
red WR bump may be a consequence of a too low S/N 



ratio. In particular, they should definitively be observed 
in NGC 5253 given the quality and depth of our spec- 
tra. Indeed, the very broad blue WR bump in region C 
of NGC 5253 strongly suggests the broad C IV A4658 emis- 
sion line, with the nebular [Fein] A4658 line on top of it 
(see Fig. [I]). We will use this line to estimate the number 
of WCE stars in this particular zone of NGC 5253. 

The line Cm A5696, much fainter than the Civ A5808 
line, is not detected in any case. Hence the emission is likely 
due to early- type WC stars (WCE). 

3. Fitting of the WR features 

As BKD08 pointed out, the fitting of the WR features, spe- 
cially in the blue bump, should be done carefully and not 
only considering one single Gaussian fit. Furthermore, the 
nebular emission lines within the bump should be prop- 
erly removed and not included in the flux of the broad 
stellar lines. Following the procedure described in BKD08, 
we performed a detailed analysis of each spectrum, check- 
ing even those spectra without a clear WR detection. We 
fitted a broad and a narrow Gaussian for the stellar and 
nebular He II A4686 lines in the blue bump and a narrow 
Gaussian for the nebular [N n] A5755 emission line ontop 
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Fig. 3. Intensity of the nebular Hell A4686 line vs. the ++ /0 + ratio (left) and the oxygen abundance (right). In both figures, 
red squares indicate objects in which the nebular Hell line but not blue WR bump are detected. Some objects have been labeled. 
The green star indicates the position of NGC 5253, taking an average value of all analyzed regions. 



a broad Gaussian for the stellar C iv A5808 line in the red 
bump. We considered the typical FWHM of the nebular 
lines to constrain the width of the nebular He II A4686 line 
and a maximum of FWHM~47 A (er=20 A) in the blue 
and FWHM~59 A (cr=25 A) in the red for the width of 
the broad stellar lines. This corresponds to FWHM^3000 
kms" 1 -1280 kms- 1 ), which is an adequate upper limit 
to the width of individual WR features (BKD08). For the 
broad Hen A4686 line, all the fits yield a <10 A except 
Tol 9 (ct=13 A). The fits to the broad C iv A5808 line were 
set to the maximum value for Mkn 1199 and Tol 9. 

Some examples of this fitting are shown in Figs. [T] (blue 
WR bump) and [H (red WR bump). The blue line traces our 
spectra, the continuum level is indicated by a dashed black 
line, and the best fit (usually composed by a broad and a 
nebular He II A4686 line in the blue bump) is shown with a 
black continuous line. These figures also show the residual 
spectrum after subtracting the best-fit model to our ob- 
served spectrum. As can be seen, the residual spectra are 
consistent with Gaussian noise plus the nebular emission 
lines. We did not consider any fits to the nebular [Fein] 
A4658 and [Anv] AA4711,4740 lines in the blue range and 
the nebular Hei A5875 line in the red range because they 
are clearly seen in the residual spectrum and, considering 
our spectral resolution, their contribution (if any) to the 
flux of the WR bumps should be negligible. We want to 
remark that this careful fitting of the broad and narrow 
lines was only performed before by BKD08, who used SDSS 
data; the majority of our spectra have both a higher spec- 
tral resolution and S/N ratio than those presented by these 
authors (compare their Fig. 2 with our Fig. [TJ. Although 
BKD08 have a more extensive galaxy sample, the additional 
strength of our study is the detailed analysis we performed 
for each region and galaxy (see Papers I and II), which 
yields a better determination of their properties (chemi- 
cal abundances, stellar population, star formation rate, col- 
ors, masses, ...). In particular, we used the bright doublet 
[O n] AA3726,3729 to estimate the total oxygen abundance, 
and not the much fainter [On] AA7319,7330 lines, which 
usually are the only [O il] lines available in SDSS data. 



Furthermore, as we will explain in Sect. 7, our galaxy sam- 
ple is more homogeneous than that used by BKD08, al- 
lowing us to perform additional assumptions about their 
star-formation history. 

As seen in Figs. [1] and [21 some spectra reveal very faint 
WR features. This is the case for region A in SBS 1319+579, 
for which a very faint blue WR bump and a nebular 
He II emission line were apparently detected, besides their 
very low S/N ratio. For the central regions of galaxies 
IRAS 08208+2818 and Tol 9, the fit of a faint broad C iv 
A5808 line also suggests WCE stars. We will consider the 
values derived from these tentative fits as upper limits to 
the actual ones. 

The broad Nin A4640 line is observed in three ob- 
jects (HCG 31 AC, Mkn 1199 and NGC 5253 C). We 
also included this line in the fitting of the data, as Fig. [1] 
shows. This broad feature is clearly evident in region C of 
NGC 5253. Incidentally, this region also shows a detection 
of the broad C iv A4656 emission line, which was also in- 
cluded in the fit of the blue WR bump. 

Table [T] compiles the results of our best fits for the neb- 
ular and broad He II A4686 and the broad C iv A5808 lines 
in those regions where these features have been detected. 
We indicate both the equivalent width and the reddening- 
corrected I/I(Hf3) ratios of these lines, including an error 
estimation for all values. We remind that these features are 
usually weak, which makes their associated uncertainties 
considerably high in some cases. 

Although we did not observe any WR feature in the 
two galaxies with the lowest metallicity of our sample 
(SBS 1211+540 and SBS 1415+437), we estimated an up- 
per limit to the broad Hen A4686 flux. These results are 
also shown in Table [1] 

4. The nebular Hen line 

As was briefly mentioned in the introduction, there is still 
some discussion about the origin of the nebular He n A4686 
line. It is generally assumed that the hard radiation field of 
WR stars produces th is nebular line in extragalactic objects 
(ISchaererlll99"l 11998ft but some authors are skeptic about 
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Table 2. O and WR stellar populations for the galaxies analyzed here in which WR features are detected. 



SV98° CH06 

Galaxy Zone WNL O c WNL O rf WCE= O' ^MBo wr£!o Wm! 



WNL+O " ""^ WNL+O WCE+O WNL WR+O WR+O 

(g) M M [*] [*] (g,h) (e) (e) (e) (i) 



HCG 31 


AC 


312 


10798 


0.028 


427 


9646 


203 


7075 


0.042 


0.028 


0.476 


0.082 


0.05V 




F2 


4.6 


200 


0.022 


7.1 


192 




192 


0.036 








0.015 


Mkn 1087 


C 


7181 


24680 


0.225 


7994 


16509 




16509 


0.326 








0.115 


Haro 15 


C 


347 


8806 


0.038 


432 


8271 




8271 


0.050 








0.037 




A 


173 


7963 


0.021 


257 


7542 




7542 


0.033 








0.021 


Mkn 1199 


C 


960 


6063 


0.137 


975 


5985 


932 


94 


0.140 


0.909 


0.956 


0.953 


0.127 J ' 


Mkn 5 


A 


2.7 


143 


0.018 


4.1 


134 




134 


0.029 








0.022 


IRAS 08208+2816 


C 


1344 


33338 


0.039 


1715 


29611 


932: 


17828: 


0.055 


0.050: 


0.543: 


0.129: 


0.051 


IRAS 08339+6517 


#1 


658 


8444 


0.072 


782 


7819 




7819 


0.091 








0.050 


POX 4 


c 


122 


4346 


0.027 


191 


4058 


113 


3464 


0.045 


0.032 


0.590 


0.081 


0.022 


UM 420 




398 


17453 


0.022 


660 


16134 




16134 


0.039 








0.021 


SBS 0926+606 


A 


32 


1987 


0.016 


54 


1897 




1897 


0.028 








0.014 


SBS 0948+532 




230 


11460 


0.020 


358 


10921 




10921 


0.032 








0.017 


SBS 1054+365 


C 


1.1 


52 


0.020 


1.7 


48.1 




1.7 


0.034 








0.023 


SBS 1211+540 




<0.27 


>14.5 


<0.018 


<0.58 


>13 




>13 


<0.043 










SBS 1319+579 


A 


6.0 


321 


0.018 


9.2 


312 




312 


0.029 








0.011 


SBS 1415+437 


C 


<1.0 


>62.8 


<0.016 


<2.4 


>57 




>57 


<0.040 










III Zw 107 


A 


309 


5567 


0.053 


420 


5010 




5010 


0.077 








0.041 


Tol 9 


C 


244 


981 


0.199 


272 


842 


127: 


37: 


0.244 


0.776: 


0.469: 


0.916: 


0.090 


ESO 566-8 




8097 


53392 


0.132 


9574 


45975 


4768 


15839 


0.172 


0.231 


0.498 


0.475 


0.073 


NGC 5253 


A 


0.73 


35.0 


0.020 


1.03 


34.2 




34.2 


0.029 








0.012 




B 


0.74 


32.0 


0.022 


1.03 


31.3 




31.3 


0.032 








0.012 






3.2 


49.2 


0.061 


4.2 


42.8 


0.65 


37.7 


0.090 


0.017 


0.154 


0.114 


0.084^ 



a We consider the solar value for the luminosity of a single WNL star to the broad Hell A4686 emission line (SV98 method). 

b We consider a metallicity-dependent value for the luminosity of a single WNL star to the broad Hell A4686 emission line (CH06 method). 

c Total number of O stars after correcting by the ionization of the WNL stars (Eq. 1 10 ft using the SV98 method. 

d Total number of O stars after correcting by the ionization of the WNL stars (Eq. HOC using the CH06 method. 

e The number of WCE stars and the WCE/ (WCE+O), WCE/WNL and WR/(WR+0) ratios were computed assuming the 

metallicity-dependent (CH06) method. 

* Total number of O stars after correcting by the ionization of the WNL and WCE stars (Eq. 1 10 I t using the CH06 method. 
9 This coincides with the WR/ (WR+O) ratio when no emission from WCE stars is detected. 

h This considers the total number of O st ars after correcting by th e ionization of the WNL and WCE stars derived using the CH06 method. 
1 Using the empirical calibration given bv lSchaerer fc Vaccal JT998T) (Eq.[TTJ. 

J Considering not only the flux of the broad Hell A4686 flux but also the Nm A4640 and Civ A4658 fluxes. 
k Number of WCE stars computed using the Civ A4658 line. See text for details. 



this (|Garnett et al.lll99ll: lGarnettll2004l) . The nebular Hell 
A4686 line may be also produced by hard X-ray radiation 
from either massive binaries or by supernova remnants 
from previous generation of stars in the range 10-50 Myr 
(iPakull fe Angebauldll986t lvan Bever fc Vanbeverenl kOOO; 



etti 

03; 



ICervino. Mas-Hesse fc Kunthl I2002TI as some observa- 
tional studies sugg es t (|Bresolin. Kennicutt fc Garnet 
19991 iGuseva et all l200fi IStasihska fc Izotovl 1200 
Thuan fc Izotovl 20051 ) , especia ll y at t he late stages of the 
star-formation burst. iGarnettl ( 20041 ) even suggests that 
as we do not yet understand the evolutionary status of 
WR stars, it is premature to predict them from stellar 
evolution models. 

Recently, BKD08 re-examined this issue using the SDSS 
data. They found that the nebular He n A4686 line was al- 
ways detected in systems with log EW ~ 2.5 as long as the 
S/N of the spectrum permitted. In the same way, not ev- 
ery galaxy where this nebular line was detected shows WR 
features, which discrepancy is more pronounced at lower 
metallicities. As lGuseva et al.l (2000) already said, in those 
galaxies (usually very low-metallicity objects) the broad 
WR feature was very probably too weak to detect. If low- 
density stellar winds at lo w-metallicity lead to a higher flux 
of Hell ionizing photons (|Smith. Norris fc Crowtherl l2002; 
lHadfield fc Crowtherl 12007 ). then normal O stars and not 



only WR stars would be expect to contribute to the Hen 
ionizing flux. BKD08 finally concluded that although at 
metallicities higher than 20% solar the ionization of Hen 
is largely produced by WN stars, at lower metallicities it is 
most likely dominated by O stars, although there should be 
a probable contribution of WN stars which is undetectable 
because of the faintness of their optical features. 

We compared the intensity of the nebular Hell 
A4686 emission line with the oxygen abundance and the 
logO ++ /0+ ratio. Notice that for NGC 5253 we analyzed 
four independent regions that are very close one to the 
other (see Lopez-Sanchez et al. 2007), and hence we also 
indicate their average values in the diagrams (green stars) . 
In Fig. [3] (left) we see that objects with higher excitation 
tend to show a higher nebular He II A4686 flux. Regions for 
which the nebular line is observed but the WR bump is not 
(HCG 31 Fl, SBS 1415+437, Tol 1457-262 A and B and 
NGC 5253 D) are dispersed throughout the diagram, but 
all have logO ++ /0 + >0. On the other hand, the intensity 
of the nebular Hen line clearly decreases with increasing 
metallicity (Fig. [3J right). Objects in the low-metallicity 
regime [12+log(0/H)<8.2] have the highest values of the 
I (He u)/ 1(11(3) ratio, and broad WR features are some- 
times not detected in their optical spectra, in agreement 
with the result found by BKD08. Hence the nebular Hen 
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A4686 emission line may arise from ionizing O stars and 
not only from WN stars at low metallicities. An analysis 
similar to that shown here involving more galaxies with 
12+log(0/H)<8.2 and new models of winds of high-mass 
stars at low metallicities are needed to confirm this assump- 
tion. Below we will only consider the broad He n emission 
line as a clear indication of WR stars within these starburst 

gal axies. 

iGuseva et al.l (|2000h did not find any nebular He n 
A4686 line for objects with 12+log(0/H)>8.13, concluding 
that high-metallicity stellar models overpredict the number 
of photons with A <228 A responsible for the ionization of 
the He + . Our careful fitting of the blue WR bump showed 
that although faint, this nebular line is present in bursts 
with higher metallicities. 



5. Computing of the WR ratios 

We estimate the number of WR stars from the luminosity 
of the broad Hen A4686 line and the broad C iv A5808 line. 
Assuming that only WNL stars contribute to the luminosity 
of the broad He II A4686 line, we can derive the number of 
this subtype WR population applying 



0.5 



WNL = N\ 



WNL 



£ ob;i (He II A4686) 
L WNL (KellMQmy 



(1) 



On the other hand, the number of WCE stars can be com- 
puted from the luminosity of the broad C iv A5808 line: 



WCE = N WC e = 



L obs (CIVA5808) 
L WC e(CTV \5808)' 



(2) 



For solar metallicities 
one single WNL star 
A4686)~1.7xl0 36 ergs" 
contributes with Lw r.n 



it is well established that 
contributes with LwNLfiien 
L , while a single WCE star 
m n ±w ce(Civ A5808)-3.3xl0 36 ergs" 1 
(jVacca fe Contilll992l). These are the quantity value s com- 
monly used (|Guseva et al.l 120001: iFernandes et al.l l2004t 
iBuckalew et al.l l2005f) and the assumptions considered 
in the theoretical models (SV98, SB99). However, there 
is increasing evidence that the WR line luminosities 



are reduced at lower metalliciti es (|Crowther fe Hadfieldl 
2006; lHadfield fc Crowtherl[2006l ICrowth cr 20071). Indeed, 
Crowther & Hadfield (2006; hereafter CH06) showed how 
reduced line luminosities at low-metallicity follow naturally 
if WR winds are metallicity-dependent. Hence we also con- 
sidered the metallicity of the star-forming regions to get a 
more appropriate estimation of the number of WR stars. 
We assumed the broad-line WNL and WCE luminosities 
given by CH06 for solar and Zq/50 metallicites, 



LwN L (HeIIA4686, Z e ) -- 
MVNL (HeIIA4686, Z Q /50) = 

Lwce(CIVA5808, Z & ) -- 
Lwce(CIVA5808, Z q /50) 



1.6 x 10 36 ergs _1 , 
2.2 x 10 35 ergs" 1 , 
2.5 x 10 36 ergs" 1 , 
= 4.0 x 10 35 ergs _1 , 



(3) 
(4) 
(5) 
(6) 



(see their Tables 3 and 4) and performed a linear fit to 
thes e values considering 12+log(Q/H 'l=8.66 and 6.96 for 
Z & (|Asplund. Grevesse fc Sauval 2005) and Z Q /50, respec- 
tively. We then adopted for each burst the WR luminosities 
computed from these relations: 

LwN L (HeIIA4686) = (-5.430 + 0.812a;) x 10 36 ergs"\ (7) 
Lwce(CIVA5808) = (-8.198 + 1.235a;) x 10 36 ergs"\ (8) 



O 
+ 



0.1 



0.05 



0.02 



WNL using l_ WNL (X4686) = f(Z) (CH06) 
- Quadratic fit: 

Y = -1 .51 1 + 0.1 1 62 X +0.91 94 X z 

X = log ( F(Broad He II X4686) / F(Hp) ) 

Y = log [ WNL / (WNL+O) ] 
with F(Hp)=100 

r = 0.943 




0.004 



0.01 0.05 

l(Broad He II M686) / l(Hp) 



0.1 



Fig. 4. Number of WNL stars (derived assuming a metallicity- 
dependence of the WNL luminosities) vs. the reddening- 
corrected flux of the broad Hen A4686 line. A quadratic fit to 
the data is shown with a dashed red line. 



with x=12+log(0/H). Table [2] compiles the number of 
WNL stars derived for each burst following both the tradi- 
tional and the metallicity-depended methods, and the num- 
ber of WCE stars in objects where the red WR bump was 
detected following the CH06 assumption. For region C in 
NGC 5253, we assumed that 



£wce(CIVA4658) = 1.71 x L W ce(C IV A5808), 



(9) 



(SV98) to compute the number of WCE within this knot. 

The total number of O stars is derived via the H/3 lumi- 
nosity, assuming Lo7v(H/3)=4.76x 10 36 erg s" 1 for a 07V 
star, N > = L(K0)/L(OTV H/J). However, the contribu- 
tion of the WR stars and other O stars subtypes to the 
ionizing flux must be considered. This is done via the the 
n(t) =07V/0 paramete r (|Vacca fc ContJll992l:IVaccalll994l: 
iSchaerer fc Va cca 1998]), which depends on the age of the 
burst. We used SV98 models (see their Fig. 21) and the 
estimated age derived for each burst to get an estimation 
of rj (listed in Col. 5 of Table [1]), and then applied 



= iV = 



-\Total 



-N v 



-.WNL 
SO 



AT nW 

jVwcEy 



(10) 



with Ql ota \ Q^ NL , Q^ CE and Qfi™ the total, the WNL, 
the WCE and the 07V numbers of ionizing photons, respec- 
tively. Following CH06, we assumed log Qn VNL =49.4 and 
log Q^ CE =49.5, while we adopted \ogQ^=49.0 (SV98). 



Q^ otal was derived from Q^ otal = iVo'Qo ■ The total 
number of O stars estimated for each burst are compiled 
in Table [21 Notice that this table includes several val- 
ues for the total number of O stars depending on the 
assumed number of WR stars (only WNL stars follow- 
ing the standard method, only WNL stars considering the 
metallicity-dependent method, or WNL and WCE stars as- 
suming the metallicity-dependent method). Table [5] also 
lists the WNL/ (WNL+O) ratios derived for the two meth- 
ods assumed here, and the WCE/(WCE+0), WCE/WNL 
and WR/ (WR+O) ratios computed using the metallicity- 
dependent method for regions in which WCE stars are de- 
tected. The number of WR stars was computed assuming 
WR=WN+WC=WNL+WCE. For galaxies in which WCE 
stars are not detected, the total number of WR stars should 
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be the number of WNL stars. However, we must keep in 
mind that because of the faintness of the broad C iv A5808 
lines (specially at low metallicities), some WCE stars con- 
tribution may be also expected, and hence the total number 
of WR stars assuming WR=WNL is just a lower limit to 

The last column in Table shows the WR/(WR+0) 
rat io obtained when a pplyin g the empirical relation given 
by [Schaerer fe Vaccal (|l998h between the intensity of the 
blue WR bump and the WR/(WR+0) ratio, 



log 



WR =a85x]o /(BlueWRBu m p)_ alL 



WR + O 



(ii) 



We assumed that the luminosity of the blue WR bump 
is the luminosity of the broad He n A4686 line, and added 
the contribution of the N in A4640 and the C iv A4658 to 
the objects for which these broad lines are detected. As 
we see, the results provided by this empirical relation are 
systematically lower than those derived using the direct 
method (Col. 5 in Table [2]) except for the three regions in 
which other broad emission lines are detected in the blue 
WR bump. Figure H plots the derived WNL/(WNL+0) 
ratio using CH06 data vs. the intensity of the broad He II 
A4686 line. As it can be seen, the best fit to the data shown 
in Fig. @] is a quadratic fit, 



y 



-1.511 + 0.1162a; + 0.9194a; 



(12) 



where x = log[J(Heii A4686)//(H/3)] assuming 7(H/3)=100 
and y = log[ WNL / (WNL+O)]. The correlation coefficient 
of this fit is r=0.943. A similar relation between the broad 
Civ A5808 intensity and the WCE/(WCE+0) ratio seems 
to be found. However, we did not attempt to fit the data 
due to the scarcity of observational points. 

6. The number of WR stars-metallicity relation 

Following evolutionary models of massive stars and popu- 
lation synthesis models (Mas- Hesse & Kunth 1991; SV98; 
Maeder & Meynet 2005), we should expect an increase of 
the number of WR stars with increasing metallicity because 
of the decrease of the minimum mass that a massive star 
needs to reach the WR phase due to the strong metallicity - 
dependence of the ste llar winds (|Vink fe de Koterl [2005) . 
Kun th fc Schildl (IT981 were the first that noticed that the 
fraction of WR stars relative to other massive stars in- 
cr eases with increasing metallicity. The analysis performed 
bv lGuseva et al.1 (|2000h also showed a tendenc y between the 
WR/ ( WR+O) ratio and the metallicity, and lPindao et al.l 
(2002) reported this dependence for high- metallicity ob- 
jects. Recent studies indicate that this correlation is cer- 
tainly satisfied (Zhang et al. 2007; Crowther 2007; BKD08). 
Besides the uncertainties and considerable high scatter of 
their data, BKD08 indicated that models of stellar evolu- 
tion including massive binaries bv lEldridee. Izzard fc Tout! 
appear to well reproduce the relation between the 



metallicity and the number of WR stars even at low metal- 
licity 

Figure [S] plots the evolution of the relative numbers of 
massive stars, WR/(WR+0), with respect to the oxygen 
abundance for our galaxy sample. Blue points represent 
WR=WNL assuming a constant luminosity for the broad 
Hen A4686, green triangles indicate WR=WNL consider- 
ing a metallicity-dependent broad Hen A4686. Dark yel- 
low stars plot WR=WNL+WCE assuming the metallicity- 
dependent of both WNL and WCE luminosities. From this 



O 

| 0.1 



DC 



0.03 



0.01 



• WR = WNL, L mBL (14686) = 1 .7x10 erg s (SV98) 
A WR = WNL, 1^^(14686) = f(Z) (CH06) 

* WR = WNL + WCE 

♦ Mkn 1087 (Va97), 
WR = WNL 



IRAS 08208+2016 
NGC 5253 0^ ^ 
HCG 31 AC i A 
* ! ! • 




SBS 1415+437 



SBS 1211+540 



NGC 5253 A,B 



7.4 



7.6 



i.O 8.2 8.4 

12 + log (O/H) 



8.6 



9.0 



Fig. 5. WR/(WR+0) vs. oxygen abundance for our galaxy 
sample. Blue points represent WR=WNL assuming a constant 
luminosity for the broad Hen A4686 emission, green triangles 
indicate WR=WNL considering a metallicity-dependent broad 
Hen A4686. Dark yellow stars plot WR=WNL+WCE assum- 
ing the metallicity-dependent WNL and WCE luminosities. 
We include the upper limit to the WR/(WR+0) ratio esti- 
mated for the very low-metallicity galaxies SBS 1211+540 and 
SBS 1415+437. Dashed lines connect the WR=WNL and the 
WR=WNL+WCE values for the galaxies for which we have 
WCE data. These galaxies have been labelled. 



figure is evident that the number of WR stars increases 
with increasing metallicity, specially in the high-metallicity 
regime. 

We should expect some scatter in Fig. [5] because of the 
strong dependence on time of the WR star relative num- 
bers during the short WR stage of the starburst. That is 
what may be happening in NGC 5253: regions A and B 
show a low WNL/ (WNL+O) ratio, region C has a very 
intense blue WR bump and a high WNR/ (WR+O) ratio, 
but no broad WR feature is found in region D even though 
all zones are very c lose. Actually, as it was explained by 
iGuseva et al.l (|2000D . a galaxy undergoing a star- formation 
burst would describe a loop in this diagram, quickly reach 
the WR/(WR+0) maximum and after that drop down to 
slightly higher metallicities because of the chemical enrich- 
ment of the surrounding interstellar medium. We consider 
that the scatter in Fig. [5] is not so evident because all ana- 
lyzed objects are undergoing an intense and very probably 
short star- formation event with an age between 3.5 and 
5.5 Myr (see Paper I) and have a very similar nature, with 
the majority of them BCDGs. 

However, Fig. [5] may suggest a possible change of 
the slope of the relation between the number of WNL 
stars and the metallicity when considering the metallicity- 
dependent broad Hen A4686 (green triangles). For galax- 
ies with 12+log(0/H)>8.2, the WNL/(WNL+0) ra- 
tio clearly increases with metallicity, but in the range 
12+log(0/H)=7.9-8.2 it seems to become constant, with a 
value of WNL/(WNL+O)-0.03-0.04 (11 points). We can- 
not confirm that tendency, however, because of the lack of 
objects in the range 12+log(0/H)=7.6-7.9, and hence more 
data in the low metallicity regime are needed to understand 
the metallicity dependence of the WNL/(WNL+0) ratio. 
This behavior was previously noticed by BKD08 and sug- 
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gested bv iGuseva et all (|2000f ). One of the main problems 
to acquire these data is that both the WR/(WR+0) ratio 
and the duration of the WR phase decrease with decreas- 
ing metallicity, so the number of galaxies with extremely 
low metallicity -12+log(0/H)<7.8- containing a WR stel- 
lar population is expected to be small. Despite the good 
spectral resolution and S/N ratio of our WHT spectra, we 
did not detect a hint of WR stars in any of the regions 
within our two lowest metallicity galaxies (SBS 1211+540 
and SBS 1415+437; their upper limits are plotted in Fig. [5]), 
only the nebular Hen A4686 line in the brightest region 
of SBS 1415+437. However, WR stars have been detected 
in the lowest met allicity galaxies known so far: I Zw 18 , 
with Z ~Z ? /50 (jlzotov et all Il997t iLegrand et ail Il997t 



iBrown et all 12002]) and SBS 0335-052E with Z ^Z^/AO 



(llzotov fe Thuanlll999l: llzotov et al1ll999t iPapaderos et al.l 
120061) . 

Ideally, we would also like to compare the WCE and 
WR=WNL+WCE corrected values between high- and low- 
metallicity objects, but this is difficult because WCE stars 
are undetected most of the time, in particular at low 
metallicities. Considering our data plotted in Fig. \E\ the 
WCE/(WCE+0) ratio also increases with metallicity in 
galaxies with 12+log(0/H)>8.2. However, the two low- 
est metallicity objects in which we detect both WNL and 
WCE stars, POX 4 and HCG 31 AC, have very similar 
WNL/(WNL+0), WCE/(WCE+0) and WR/(WR+0) 
ratios (see Table [2]) despite their relatively different oxy- 
gen abundance, 12+log(O/H)=8.03 and 8.22, respectively. 
Again, we state that more data of low-metallicity galaxies 
with clear detections of both the blue and red WR bumps 
are needed to understand this behavior. 

Another interesting information from Fig. [5] is that for 
all galaxies with a detection of WCE stars the derived 
WR/ (WR+O) ratio is considerably higher than when only 
WNL stars are considered. The contribution of these WCE 
stars is smaller at lower metallicities, but even for POX 4 
(the lowest metallicity galaxy where both WNL and WCE 
stars are detected) the WCE stars contribute to the 37% 
to the total WR number. This as well as the analysis of 
the WCE/ WNL ratio that we will explain below strongly 
suggests that WCE stars are not only present at low metal- 
licities but also make an important contribution to the total 
number of WR stars. 



7. Comparison with the theoretical models 

The first models that tried to quantify the popula- 
tion of WR stars in sta r bursts were those presented by 
I Arnault. Kunth fe Schill (jl989t) . who found that the star- 
formation event only occurs during a short period of time 
compared with the typical lifetime of the massive stars 
and that the intensity of the blue WR bump with re- 
spect H/3 decrea ses with decreasing metallic ity. The mod- 
els presented by Mas-Hesse fe Kunthl ((1991), updated by 
ICerviho fe Mas-Hessd (|1994L were the first models that de- 
termined the number of WR stars. However, their method 
had several flaws because they did not separate the stel- 
lar emission lines within the blue WR bump, which as we 
see are also blended with the nearby Fe, He and Ar lines 
from the emission of the nebular gas. W ith the aim of solv- 
ing this problem, iKriieer et al.l (|l992l ) separately synthe- 
sized the broad Hen A4686 and Ciii/iv A4650. iMevnetl 
(|l995l) studied the effect of changing the star-formation 



rate, the IMF, the age and the metallicity in the massive 
star populations. Alth o ugh t he detailed models provided by 
iLeitherer fe Heckmanl (|1995l) included a good study of the 
O and WR populations, they failed to reproduce some ob- 
servational features. The evolut i onary models presented by 
iGarcfa- Vargas. Bressan fe Dfazl ()1995[ ) had the same prob- 
lem. 

Schaerer & Vacca (1998) (hereafter SV98) constructed 
evolutionary synthesis models for young starbursts using 
stellar evolution models, theoretical stellar spectra and a 
compilation of observed line-strengths from WR stars. They 
explicitly distinguish between several WR stars (WN, WC, 
WO) whose relative frequency strongly depends on the 
metallicity, and also consider O and Of stars independently. 
The SV98 models compute the number of O and WR stars 
produced during the starbursts and predict the intensities 
of the UV and optical emission lines of both the main neb- 
ular lines (H and He) and the broad stellar lines (He, N and 
C) as a function of several parameters related to the star- 
formation episode. These models provide the most reliable 
determinations to date. The basic framework of the SV98 
models was recently included in the STARBURST 99 stel- 
lar population synth esis model (ILeitherer et al.lll999f) with 
refined spectra fro m ISmith et al.l (l2002h. The last release 
of the SB99 code (jVazauez fe Leithererl 120051 ) is available 
on-lineQ. 

We considered both the SV98 and SB99 models to com- 
pare their predictions with our observational data. In both 
cases we assumed an instantaneous burst with a Salpeter 
IMF slope a=2.35 and M up =120 M Q with three different 
metallicities, Z=0.004, 0.008 and 0.020 (Z Q ), which corre- 
spond to an oxygen abundance -in units of 12+log(0/H)- 
of 7.96, 8.26 and 8.66, respectively. The SV98 models con- 
sider Mi ow =0.8 Mq and SB99 models assume Mi ow =0.1 
Mq. Although this does not affect to the predictions for 
the massive stars, the SB99 models were created assuming 
Padova tracks in which asymptotic giant branch (AGB) 
phases with the same met allicity of the gas are included 
(jVazauez fe Leithererll2005r) . The total stellar mass created 
in the starburst (a normalization factor) was set to 1 Mq 
and 10 6 M Q for SV98 and SB99 models, respectively. 

Besides that the SV98 models are the most reliable 
WR galaxy models to date, we remind the reader that 
nowadays there is a better understanding of the properties 
of the WR stars than when the SV98 models were released. 
The main refinement is the metallicity-dependence of 
the Hen A4686 l uminosity (CH06) and t he influence 
of the rotation (jMevnet fe Maeden 120051) and wind 
loss and binary evolution (Ivan Bever fe Vanbeverenl 



2003; 



Vanbeveren. van Bever fe Belkusl 



Eldridge. Izzard fe Tout! 120081) in the stellar 



2007; 
tracks 

followed by massive stars. Stellar rotation is expected to 
predict a longer WR stage, lower WR star masses and 
thus larger WR populations. BKD08 considered some 
variations in the WR lines fluxes with the metallicity 
to update the SV98 models, but these models were not 
published anywhere. These authors also remarked the 
importance of the star-formation histories of the galaxies 
hosting the WR-rich starclusters, combining their models 
with the predictions giv en by both the SB99 and the 
iBruzual fe Charloll (|2003l) codes to generate a wide range 
of models that depend on many parameters. New models 



http:/ /www.stsci.edu/s 



e/starburst99/ 
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-W(Hp) [A] -W(Hp) [A] 

Fig. 6. Intensity (left) and equivalent width (right) of the nebular Hell A4686 line vs. W(Rl3) and their comparison with the 
predictions given by SV98 models. Three different metallicities, Z=0.004, 0.008 and 0.020 (Z@) are plotted with dashed red, 
continuous blue, and dashed-dotted green lines, respectively. Data colors (red, blue and green) indicate the model that best 
matches the metallicity of the object. 



of evolving starbursts showing WR features that include 
all these aspects are absolutely needed. 

Figures IBTfTUl compare the predictions of the SV98 
and SB99 theoretical models with our observational data. 
As the luminosity of the WR features and their relative 
strength compared to other spectral properties depends 
on age, metallicity, star-formation history, and the IMF, 
the comparison of the observational data with the mod- 
els should be carefully done. In all these figures the green 
dashed-dotted line represents the model predictions for so- 
lar metallicities (Z=0.020), the continuous blue line is the 
prediction of the Z=0.008 model and the red dashed line 
indicates the predictions of the Z=0.004 models. The color 
of the observational points is chosen according to the model 
that better matches their metallicity. We always represent 
in the a;-axis the H/3 equivalent width, W(Hf3), which is a 
good indicator for the age of the most recent massive star- 
formation event. Indeed, the hydrogen- ionizing flux of a star 
cluster gradually decreases as the most massive stars disap- 
pear with time, and hence the absolute value of the width 
of H/3 decreases with time (see Papers I and II). BKD08 
questioned the validity of this parameter for a large sample 
of WR galaxies spanning a wide range in metallicities and 
sizes, moreover if the spectrum covers a big region in the 
galaxy including young and old stellar populations. This 
is not typically our case, because all the spectra analyzed 
here were mainly obtained in brightest Ha-selected regions 
within the galaxies (see Paper II) in which the underlying 
old stellar population is not dominating the optical light. 
Hence, while the spectroscopic SDSS data used by BKD08 
randomly selected a very heterogeneous sample of objects 
with a fixed fiber size (3"), we first carefully analyzed the 
position of the brightest H II regions within our sample of 
WR galaxies (Paper I) and then extracted the spectroscopic 
data using efficient apertures (Paper II). Furthermore, both 
the analysis of the Ha imagery and the spectra indicate 
that all objects have experienced an almost instantaneous 
star-formation event that happened between 3 and 6 Myr 
ago. Although some few cases of BCDGs seem to indicate 
continuous star formation for some few Myr (i.e., Mkn 209, 
see Perez-Montero & Diaz, 2007), the majority of the anal- 



yses of similar low-metallicity galaxies [i.e., Guseva et al. 
2000; Cairos et al. 2001, 2007; Gil de Paz & Madore 2005] 
support the instantaneous star-formation scenario. Hence 
we consider the comparison of the observational data with 
a model of a single star-formation event appropriate for our 
analysis, and that W(Rf3) is a reliable indicator of the age 
of the starbursts. 

Figure H2 compares the predictions of the nebular He II 
A4686 following the SV98 models with our data. The di- 
agrams plot the intensity (left panel) and the equivalent 
width (right panel) of this nebular line as a function of the 
H/3 equivalent width. Notice that several high-metallicity 
objects -green dots- with W(Hft)> —60 A are displaced 
somewhat to the right of the model predictions in both di- 
agrams. We think that this effect is a consequence of a prob- 
able underestimation of the W (H/3) in these objects, which 
have more complex star-formation histories and a consid- 
erable contribution of the old underlying population than 
the low- and intermediate-metallicity objects. We should 
keep this in mind when comparing observed and predicted 
properties for high-metallicity objects in the next figures. 
In particular, we should expect a much more complex star- 
formation history in a spiral galaxy than in a BCDG. In our 
study, only Mkn 1199, Mkn 1087, ESO 566-8, ESO 566-7, 
IRAS 08339+6517 and Tol 9, all with 12+log(0/H)>8.4, 
show clear complex star-formation histories, although to- 
day all of them are hosting a strong starburst. These objects 
are plotted with green dots in Figs. IBTTLUl Perhaps their 
W(R(3) have been slightly underestimated, but the com- 
parison of their observed emission line properties with the 
pre dictions given by the photoionization models provided 
by IStasihska. Schaerer fc Leithererl (|200ll) (see Fig. 10 in 
Paper IV) indicates that the underestimation of the W(H(3) 
value is not very important (10-20 A at almost). Hence, 
although the W (H/3) derived in high- metallicities objects 
may be somewhat underestimated, we think that it is not 
the case for low- and intermediate-metallicity objects be- 
cause (i) the contribution of the underlying stellar popula- 
tion, although it does exist, is not dominating in the strong 
star-formation bursts (see Sect. 4 in Paper I and its Figs. 37- 
39) and (ii) the ages of the star-forming bursts obtained 
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Fig. 7. Comparison of the observed intensities (upper row) and equivalent width (lower row) of the broad He n A4686 emission 
line with the predictions given by SV98 (left) and SB99 (right) models. Symbols are the same as in Fig. [5] 



from W(Hf3) agree with those determined using indepen- 
dent methods (see Sect 4 in Paper II and its Fig. 39). 

From Fig. [6] (left) is also evident that galaxies with 
higher mctallicities show lower Hell A4686 intensities. As 
we already discussed, this line is not or only very weakly 
detected in our highest metallicities objects. Both the He II 
A4686 intensity and its equivalent width (Fig. |5l right) 
decrease with increasing W^H/?). The model predictions 
by SV98 generally agree well with our data, but perhaps 
they overpredict the intensity of th e nebu lar He n line at 
high metallicities, as iGuseva et all (|2000D previously no- 
ticed. These authors also pointed out that objects with de- 
tected and nondetected WR features mingle indistinctly in 
these plots, and hence not only WR stars but also O stars 
may be contributing to the ionization of the He + . Precisely 
this was one of the conclusions reached by BKD08, as we 
already discussed in Sect. 4. 

Figure [7] compares the observed intensities and equiva- 
lent width of the broad He n A4686 emission line with the 
predictions given by SV98 and SB99 models. Again the the- 
oretical predictions agree well with the observational data, 
but we observe that the solar-metallicities model clearly 
overpredict the broad He II A4686 intensities and equivalent 
widths by a factor 4-8. Although part of this disagreement 
maybe a consequence of the more complex star-formation 



histories of the high-metallicity galaxies (and hence an in- 
crease of the IF(H/3) because of the higher importance of 
older stellar populations), we consider that some of the 
disagre ement between pred ictions and observational data 
is real. iPindao et al.l (|2002D also reported that theoretical 
models overestimate the WR fluxes in solar and super-solar 
metallicity objects. 

Actually, we hold that the agreement we see at low 
metallicities may be just incidental. The SV98 and SB99 
models consider the solar values of the intensities of the 
broad WR lines independently of the metallicity, but as we 
explained before (Sect. 5) they are drastically reduced at 
low metallicities. Therefore, the predictions given by these 
models should not agree with the observational data at 
low metallicities. Hence, models using the same parame- 
ters as SV98 and SB99 do agree but including metallicity- 
dependent WR line luminosities would underpredict the 
observed broad Hen A4686 intensity (and its equivalent 
width). For example, considering the formulae presented 
before (Eq. O, the model with 12+log(0/H)=7.96 (red 
dashed line in the figures) would predict broad He n in- 
tensities that are 65% of the intensities shown in Fig. 
which makes their equivalent widths also much lower than 
those seen in the figure. 
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line with the predictions given by SV98 (left) and SB99 (right) models. Symbols are the same as in Fig. [B] 



The same situation is found when comparing the ob- 
served broad C iv A5808 emission line with the model pre- 
dictions (Fig. |5]). For example, for the central burst in 
POX 4, the predicted /(Civ A5808)//(H^) is -0.015 fol- 
lowing the SV98 models, quite in agreement with the actual 
value observed in this object (0.014). Considering the effect 
of lower WR luminosities at lower metallicities (Eq. [SJ , the 
value predicted by the models should be ~0.01, which is 
around 70% of the observed value. Hence, the SV98 and 
SB99 theoretical models are indeed underpredicting the in- 
tensities of the broad WR lines at low metallicities. This 
result agrees with that reported by BKD08 and previously 
noticed by other authors (e.g. Legrand et al. 1997; Guseva 
et al. 2000). The best explanation of this discrepancy seems 
to be that the binary channel for producing WR stars is im- 
portant at low metallicities, so that single star models are 
not efficient enough to reproduce the observed values. The 
inclusi on of binaries in the models led t o a prolongated WR 
phase (jvan Bever &; Vanbeverenl l2003fl and to an increase 
in the WR/(WR+Q1 ratio (lEldridee et al.l[2008l) . The bi- 
nary channel ingredient should be taken into account in new 
models of evolving starbursts showing WR features, as it 
has not bee n included before, althou gh there is some work 
m process (|Eldridee fc Stanwavl [20091) . These new models 



should also consider the metallicity dependence of the WR 
luminosities. 

Besides this problem, from Fig. [S] is also clear that the 
predictions given by the SB99 models are worse than those 
predicted by the SV98 models. The SB99 predictions quite 
disagree for low- and intermediate-metallicity objects. Even 
without considering the diminution of the WR luminosi- 
ties, the SB99 models predict /(Civ A5808)/7(FL3)~0.2 for 
POX 4, which is only the 14% of the observed value. Hence 
we suggest that we detect more WCE stars than those pre- 
dicted by the models, which difference is even more pro- 
nounced at lower metallicities. 

Let us now compare the observed WR/(WR+0) ra- 
tios with the theoretical predictions (Kig. |9j) . The numbers 
of WNL and WCE stars predicted by the models do not 
depend on the luminosities of the broad emission features 
assumed for a single WR star (only the predicted WR 
luminosities depend on them) and hence the comparison 
with our observational data should be still valid. The SB99 
models better reproduce the WNL/(WNL+0) ratio than 
the SV98 models (Fig. EE top row): the WNL/(WNL+0) 
ratio predicted by the SB99 models agrees well with the 
observed values for all metallicities, while the SV98 mod- 
els underestimate the W(Hf3) in which the maximum of 
the WR/(WR+0) ratio is expected, which is very evident 
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at low metallicities (differences >100 A). If our observa- 
tional points were affected by underlying absorption, this 
effect would be even higher. However, both kinds of models 
strongly disagree with the observations at high metallicities 
when compared with the observed WCE/(WCE+0) ratio 
(Fig. [SI central row). The disagreement is particularly im- 
portant (more than one order of magnitude) for Mkn 1199. 



The lower row in Fig. shows the comparison of the 
derived WR/(WR+0) ratio with the theoretical predic- 
tions given by SV98 and SB99. In both cases we con- 
sidered the ratios involving WNL and WR=WNL+WCE 
stars. The number of WCE stars estimated for Tol 9 
and IRAS 08208+2816 is tentative, so their derived 
WR/(WR+0) ratio is an upper limit to the real one. As 
we commented before, there is a disagreement between the 
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Fig. 10. Comparison of the derived WCE/WNL ratios with the predictions given by the SV98 (left) and SB99 (right) models. 
Symbols are the same as in Fig. [6] 



SV98 predictions and the observations at low metallicities, 
as these tend to show higher iy(H/3) values for the de- 
rived WR/(WR+0) ratio. If WCE are typically present 
in those bursts, the models would also underpredict the 
WR/(WR+0) ratio. Considering again the low-metallicity 
galaxy POX 4, the SV98 predictions for its WR/(WR+0) 
ratio are ~4 times smaller than the observed value consid- 
ering the detection of both WNL and WCE stars in this 
object. This difference is also observed in objects with in- 
termediate metallicites, although SV98 models work better 
in the solar-metallicity range. 

On the other hand, the SB99 models seem to overpre- 
dict the WR/ (WR+O) ratio at low luminosities. As we saw 
that the observed WNL/(WNL+0) ratio agrees well with 
the predictions in SB99 models, we may interpret this result 
in two ways: (i) models are overpredicting the number of 
WCE stars or (ii) we are not detecting WCE stars in these 
low-metallicity objects. In the only low-metallicity galaxy 
where we detect both WNL and WCE stars the agreement 
is not bad: the SB99 models account for around 80% of the 
WR population detected in POX 4. As we see, in all regions 
in which both subtypes of WR stars are detected the mod- 
els underestimate the total number of WR stars. Hence the 
problem seems to be a consequence of the predictions of the 
number of WCE stars, but no the number of WNL stars. As 
we will see below, the derived WNL/ WCE ratio computed 
for these objec ts is higher than that predicted by more 
recen t models (jMevnet k Maederl 120051: lEldridee k Vinkl 
|2006|) and also higher than that empirically found counting 
individual WR stars in nearby galaxies. 

8. The WCE/WNL ratio 

Table [2] compiles the WCE/WNL ratio computed in the 
galaxies for which both WNL and WCE stars have been 
detected. Figure ITU1 compares the WCE/WNL ratios with 
the results given by the SV98 and SB99 models. The SV98 
models clearly predict an increase of the maximum of the 
WCE/WNL ratio with increasing metallicity, but that is 
not seen in the SB99 models. Actually, the predictions given 
by these two kinds of models clearly disagree. The SV98 
models seem to better reproduce the observed ratios at 
low and intermediate metallicities, but because of all the 



caveats about the predictions of the numbers of WNL and 
WCE stars we consider that the SV98 and SB99 models are 
not valid. 

A much more interesting analysis is the metallicity- 
dependence of the WCE/WNL ratio. Figure Qj] shows the 
WCE/WNL ratio as a function of the oxygen abundance 
for our starbursting galaxies (blue circles) and its compari- 
son with the values obtained for nearby spiral (red squares) 
and irregular (yellow stars) gal axies for which t he WR con- 
tent has been studied in detail (|Crowtherll2007D . Except for 
the par ticular case of the starbursting Local Group ga laxy 
IC 10 (|Massev k Holmes! [200l ICrowther et Zl 120031) . for 
which the true WR populations are still somewhat uncer- 
tain, the WCE/WNL ratio clearly increases with the metal- 
licity. That trend also agrees with the SV98 and SB99 
models predictions (not shown in Fig. [11]), for example, 
the maximum of the WR/(WR+0) ratio in SV98 mod- 
els increases from 0.025 to 1 between metallicities of 7.36 
and 8.96 - i n unit s of 12+log(0/H)-, respectively. Following 
ICrowtherl (|2007T ). we plot in Fi g. ITT1 the predictions give n 
by the most recent WR models. [Mevnet k Maederl ()2005l) . 
which allow for rotatio nal mixing but without a WR metal- 
licity dependence, and lEldridge k Vinkl (|2006[ ). which ne- 
glect rotational mixing but include a metallicity-scaling 
for W R stars. At high metallicities, the Eldrid ge k Vinkl 
(2006) models agree well with previous observations and 
also with our value for Mkn 1199, a spiral galaxy experienc- 
ing a strong star-formation event because of the interaction 
with a nearby companion (see Papers I and II). 

However, the situation is quite different in the 
intermediate- and low-metallicity regime. Previous data of 
spiral and irregular galaxies agree well with both kind of 
models, but except for region C within NGC 5253, they 
clearly underestimate the WCE/WNL ratio in our sam- 
ple of starburst ga laxies. That is the same result that 
iGuseva et al.l (|2000l) found in their analysis. They assumed 
that the WC / WN ratio empirically derived by counting in- 
dividual WR stars in nearby galaxies corresponds to con- 
tinuous star formation, that should be ~0 for galaxies with 
12+log(0/H)<8.1. But low-metallicity galaxies are expe- 
riencing an instantaneous starburst, and hence their de- 
rived WCE/WNL may be as high as -3 for Z=0.004 - 
12+log(0/H)~7.96- according to SV98 models. This ob- 
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Fig. 11. Comparison between the observed WCE/WNL ratio 
and the metallicity for the starburst galaxies analyzed in this 
work (blue circles). We include the WC/WN ratios derived for 
nearby spiral (red squares) and irregular (yellow stars) galaxies 
(Massey & Johnson, 1998; Crowther et al. 2003, 2007; Schild et 
al. 2003; Hadfield et al. 2005; Hadfield & Crowther 2007). Several 
regions within M 33 are shown because of the strong metallicity 
gradient in this galaxy. The evolutionary model predictions by 
Meynet & Maeder (2005, black continuous line) including rota- 
tion and Eldridge & Vink (2006, green dashed line) including a 
metallicity- dependence of the WR winds are also shown. 



servation also agrees very well with the high ratio de- 
rived for IC 10, WCE/WNL~1.25. We therefore conclude 
that the differences in the WCE/WNL ratio found between 
starbursting low-metallicity dwarf galaxies and the typi- 
cal nearby star-forming galaxies are a consequence of their 
very different star-formation histories, and that new models 
should definitively include this component. 

However, we must keep in mind the effect of the slit 
position and the size of the aperture on the derived WR 
properties, which is particularly evident in nearby objects 
such as IC 10 or NGC 5253. Indeed, we see that different 
regions within these galaxies have very different WR pop- 
ulations. To investigate this issue, we plot in Fig. [12] the 
metallicity-dependence WNL/(WNL+0) ratio as a func- 
tion of the distance to the galaxies (see Table 1 in Paper I) . 
The distance range of our sample galaxy is between 4 Mpc 
(NGC 5253) and 237 Mpc (UM 420), with the majority of 
the objects within a distance of 100 Mpc. We do not see 
any clear correlation in this diagram, so we can discard any 
distance effect in our conclusions. We remind that we did 
not consider all flux of the galaxy with a fixed aperture, but 
analyzed the spectra of the brightest H II regions found in 
each galaxy using our Ha maps (Paper I) and choosing the 
size of the aperture according to these regions (Paper II). 
Hence the effects of dilution of the WR features (includ- 
ing areas in which the WR stars are not presented) were 
minimized in our analysis. 



9. The WR dependence on the IMF 

The detection of WR stars of both WN and WC subtypes 
and the derived WC/WN ratio provides strong constraints 
for stellar evolution models through the determination of 
the upper cut-off limit of the IMF. As WC stars are diffi- 
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Fig. 12. WNL/(WNL+0) ratio vs. the distance to the galaxies. 
Notice that we also indicate the region within the galaxy in 
which the WR features are detected. The average value of the 
WNL/(WNL+0) ratio for the BCDG NGC 5253 is plotted with 
a yellow star. 



cult to observe in low-metallicity environments, these stud- 
ies are usually performed in metal-rich objects. Actually, 
the slope of the IMF and its dependence on metallicity is 
still a controversial issue. Som e studies support a Salpeter- 
like IMF at solar m e tallicities (|Bresolin fc Kennicuttll 2002t 
Pindao et all l2002t iFernandes et al.l l2004t iBresolin et all 
20051 ). but other observations hav e questioned t h is as- 
sum ption at different m etallicities. iSchaerer et al.l fl999) 
and iHuang et al. (I1999D sugg e sted a flatter IMF at lower 
metallicities. Fcrnande s et al. (|2004l ) underlined a depen- 
dence of the IMF slope on galaxy metallicity, in which low 
metallicity galaxies show a Salpeter-like IMF while high- 
metallicities galaxies either show a steeper IMF or expe- 
rience an extende d burst. A s o mewh at similar result was 
recently found by Zhang ct al. (2007), who found that the 
IMF slope increases with increasing metallicity, with the 
slope index ranging from a ~1.0 for Z=0.001 to a ~3.3 for 
Z=0.020. 

BKD08 pointed out that taking into account the con- 
siderable number of parameters that should be considered 
when modeling starbursts with significant WR populations 
(age of the burst, star-formation history, metallicity, de- 
pendence of the WR features on stellar rotation and wind 
loss, inclusion of the binary channel...) it is not a good idea 
to include the slope of the IMF as an additional parame- 
ter. Just as an exercise, we compared the predictions given 
by SV98 models for an instantaneous burst with the same 
properties we explained before but changing the slope of 
the IMF, a=1.00, 2.35 (Salpeter) and 3.00. Figure QJ plots 
the observed WR/(WR+0) ratio versus the oxygen abun- 
dance and the maximum of the WR/(WR+0) predicted by 
these three models. The observational data plotted in this 
figure are the same as in Fig. [5] the WNL and WCE stars 
numbers were computed assuming a metallicity-dependent 
luminosity of the WR features, as we explained before. As 
we see, models with flatter IMFs tend to predict higher WR 
values and at all metallicities, so that a flat IMF implies a 
relatively large population of massive stars, leading to a 
larger WR population. 
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Fig. 13. WR/(WR+0) vs. metallicity and the predictions of 
SV98 models assuming three different IMFs, a=1.00 (flatter, red 
dashed line), 2.35 (Salpeter, dashed-dotted purple line) and 3.00 
(steeper, blue continuous line). As in Fig. [5] green triangles indi- 
cate WR= WNL considering a metallicity-dependent broad He II 
A4686, and dark yellow stars plot WR=WNL+WCE assuming 
the metallicity-dependence of both WNL and WCE luminosities. 
Dashed lines connect the WR=WNL and the WR=WNL+WCE 
values for the galaxies for which we have WCE data, these 
objects have been labelled. We include the upper limit to 
the WR/(WR+0) ratio estimated for the very low-metallicity 
galaxies SBS 1211+540 and SBS 1415+437. 



From Fig.[T3]is evident that our data do not support the 
hypothesis that low-metallicity objects have natter IMFs. 
Actually it is quite the opposite, the model with a=3.00 
seems to agree better with the data. The assumption of 
metallicity-independent WR luminosities would move our 
data points to lower WR/(WR+0) values, and hence fa- 
vore an even steeper slope for the IMF at low metallic- 
ities. However, we notice that the majority of the data 
points assumes WR=WNL, but we have seen that the WCE 
stars should have a non-negligible contribution to the total 
WR star population even at low metallicities. Actually, the 
model with a Salpeter-like IMF agrees very well with the 
observed WR/(WR+0) value of the objects for which we 
have detected both WNL and WCE stars. We then con- 
clude that the assumption of a Salpeter-like IMF is valid 
and that the metallicity-dependence on the WR luminosi- 
ties is probably playing a key fact or in the discrepan cies 
between data and models found by IZhang et al.l (|2007l ). 

A final rem ark about IMFs in star-forming galaxies. 
Recent studies ()Weidner fc Kroupal 120051 I2006D have re- 
vealed that the total stellar population of the new formed 
stars of all young star clusters within the same galaxy 
follows a distribution function steeper than the canoni- 
cal IMF in the high-mass regime. The so-called integrated 
galactic initial mass function (IGIMF) deviates increas- 
ingly from the underlying canonical IMF with decreasing 
star-formation and total galaxy mass. Blue Compact Dwarf 
galaxies, such as the typical low-metallicity objects in which 
WR stars are detected, should have a very different IGIMF 
than spiral galaxies like the Milky Way. This parameter 
should be therefore also considered when modelling the 
evolution of the stellar populations in starbursts and star- 
forming galaxies. 



10. Conclusions 

We presented a detailed analysis of the broad stellar fea- 
tures originated by winds of Wolf-Rayet stars in a sample of 
starbursting galaxies previously classified as WR galaxies. 
The photometric properties of these galaxies were compiled 
in Paper I, while their spectroscopic analysis was shown in 
Paper II. The metallicity of these galaxies lies between 7.58 
and 8.75 -in units of 12+log(0/H)-. Our goals are to lo- 
cate the WR-rich stellar clusters, derive the number of O, 
WN and WC stars and compare the results with previous 
observations and the predictions of evolutionary synthesis 
models. Our main results are the following: 

1. The blue WR bump, mainly composed by the broad, 
stellar He n A4686 emission line and attributed to WN 
stars, is unambiguously detected in 12 regions and it is 
probably observed in eight other regions. Aperture ef- 
fects and/or the position of the slit within the starburst 
are clearly playing a fundamental role in the detection 
of the WR features. We consider that multi-aperture 
spectroscopy or narrow-band imagery using filters sen- 
sitive to WR features and the adjacent continuum later 
followed by slit spectroscopy may be more appropri- 
ate for this kind of analysis. The advent of the optical 
3D spectroscopy is already showing the localization of 
WR features with star-forming regions (i.e. Kehrig et al. 
2008; Garcfa-Lorenzo et al. 2008; James et al. 2009,2010; 
Lopez-Sanchez et al. 2010). 

2. The broad C IV A5808 emission line (red WR bump), at- 
tributed to WC stars, is clearly detected in four regions, 
and it is probably observed in two other objects. 

3. We have performed a detailed fitting of the spectra con- 
sidering the broad stellar and narrow nebular emission 
lines in both the blue and the red WR bumps. In the 
blue WR bump we usually found a contribution of the 
nebular and broad Hen A4686 components. 

4. The nebular He n A4686 is detected in 17 regions and it 
is probably observed in six other knots. We confirm the 
trend that galaxies with lower metallicities have a higher 
/(nebular He II A4686) / 1(11/3) ratio. Some of our lowest- 
metallicity regions show the nebular but not the broad 
Hell emission line. The SV98 models seem to under- 
predict the intensity of this emission line in low metal- 
licity objects. Both findings agree with the conclusions 
reached by BKD08, who suggested that O stars and not 
only WR stars contribute to the ionization of the He + 
at low metallicities. 

5. We used the intensity ratios of the broad Hen A4686 
and Civ A5808 lines to derive the number of O, WN 
and WC stars within each region. We have considered 
two methods: (i) assuming constant luminosities for the 
WR features, as it is commonly done and it used in theo- 
retical models, and (ii) assuming metallicity-dependent 
luminosities for the WR features, as both recent ob- 
servations and new WR stars models indicate (CH06). 
This second method gives higher number of WR stars 
at lower metallicities. We derive an empirical estima- 
tion of the WNL/(WNL+0) ratio using the intensity 
of the broad Hen A4686 line assuming the metallicity- 
dependent of the WR luminosities. 

6. As expected following theoretical models, the total num- 
ber of WR stars increases with increasing metallic- 
ity. However, we observe a possible flattening of the 
WR/(WR+0) ratio for objects with 12+log(0/H)<8.2, 
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because all the low-metallicity galaxies show a constant 
value of ^0.03 - 0.04, but more data are needed to con- 
firm this trend. We also conclude that the contribution 
of the WCE stars is not negligible at low metallicities. 

7. The comparison of the observational data with the the- 
oretical predictions given by the most accurate models 
available to date, SV98 and SB99, is difficult because 
of all the parameters involved, namely, age, metallic- 
ity, star-formation history, assumption of the IMF and 
WR stars properties (variation of the WR luminosities 
with the metallicity, effect of star rotation, contribution 
of the WR binnary channel). New models including all 
these factors are absolutely needed to perform an ap- 
propriate comparison with the observational data. The 
available SV98 and SB99 models seem to work better at 
higher metallicities, but more complex star-formation 
histories arc needed to explain the observed WR lu- 
minosities and ratios in this regime. At low metallici- 
ties the SV98 and SB99 models fail to reproduce the 
observed WR intensities because these models do not 
consider the decreasing of the WR luminosities with 
decreasing metallicity. Our data agree with the SV98 
models considering a Salpeter-like IMF for all metal- 
licitcs. 

8. Finally, the WCE/WNL ratios observed in our star- 
bursting galaxies are quite different than those empiri- 
cally found in nearby star-forming and irregular galax- 
ies. We consider that this effect is a consequence of the 
very different star-formation histories that these objects 
are experiencing. 
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